The survival probability of a branching random walk in presence of an absorbing wall 
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A branching random walk in presence of an absorbing wall moving at a constant velocity v 
undergoes a phase transition as v varies. The problem can be analyzed using the properties of the 
Fisher-Kolmogorov-Petrovsky-Piscounov (F-KPP) equation. We find that the survival probability of 
the branching random walk vanishes at a critical velocity v c of the wall with an essential singularity 
and we characterize the divergences of the relaxation times for v < v c and v > v c . At v — v c 
the survival probability decays like a stretched exponential. Using the F-KPP equation, one can 
also calculate the distribution of the population size at time t conditionned by the survival of one 
individual at a later time T > t. Our numerical results indicate that the size of the population 
diverges like the exponential of (v c — v)~ 1 ^ 2 in the quasi-stationary regime below v c . Moreover for 
v > v c , our data indicate that there is no quasi-stationary regime. 



There has been a long-standing interest in problems 
of non-equilibrium critical phenomena and phase tran- 
sition into absorbing states [H H (directed percolation 
or contact processes [E HJ, reaction diffusion problems 
[E @1)- The goal of the present letter is to analyse one 
such problem, which has been already considered by sev- 
eral authors in the mathematical literature [1, [E EE] : 
a branching random walk in presence of an absorbing 
moving wall at a constant velocity v. For such a prob- 
lem, there is an obvious absorbing state (all the particles 
are absorbed by the wall) and in this letter we analyse 
the critical behaviour of the survival probability as the 
velocity v of the wall varies. Our approach is based on an 
analysis of the Fisher-Kolmogorov-Petrovsky-Piscounov 
(F-KPP) equation [H El El EI EE 03 El which is 
known to be related to branching random walks [IE EH • 
We also show that the F-KPP equation allows us to cal- 
culate the properties of the quasi-stationary regime. 

In absence of an absorbing wall, a one dimensional 
branching random walk has a number of descendants 
which grows exponentially with time and these descen- 
dants occupy a region which spreads [IE Ell m space 
at a known velocity v c (which can be calculated from 
the knowledge of the branching and hopping rates of the 
walk). In presence of an absorbing wall moving at a 
constant velocity v, there is a competition between this 
exponential growth of the number of descendants and 
the absorption by the wall. One quantity of interest is 
the survival probability Q(x,t), for a particle starting at 
t = Q at distance x from the wall, to have at least one 
surviving descendant at time t. If v > v c (i.e. if the 
wall moves faster than the spreading velocity of the pop- 
ulation in absence of the wall), Q{x,t) — > as t — ► 00. 
In this letter we study the phase transition which takes 
place at v = v c when v varies Q. We will show that for 
v < v c , the survival probability Q(x,t) has a non-zero 
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FIG. 1: Branching random walk in presence of an absorbing 
wall. The circles represent branching events and the crosses 
absorption by the wall. 



long time limit which vanishes with an essential singular- 
ity as v — * v c . We will also show that the characteristic 
time scale r (on which Q(x,t) converges exponentially 
to Q{x,oo)) diverges like \v — v c \~z below v c and like 
\v — Vc}^ 1 above v c . At v = v c , Q{x,t) decays like a 
stretched exponential as already proved by Kesten Q- 

Our motivation for the problem in presence of a moving 
absorbing wall comes from the interest for simple mod- 
els of evolution [19|, |2fJ proposed recently to study the 
in vitro evolution of a population of DNA molecules se- 
lected according to their ability to bind to some target 
proteins [2l[. In these models, the position of an indi- 
vidual represents its adequacy to the environment (or 
fitness or ability to bind), meaning that when selection 
acts, the individuals with lowest x are eliminated first. 
The branching random walk has the effect of modeling 
the growth of the population and the mutations through 
diffusion along the fitness axis. In the models studied 
in [22|, HE H3 j the size of the population is kept fixed : 
selection means that whenever new individuals appear in 
the population, the ones with lowest x are eliminated in 
order to keep the size of the population constant. Here 
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the role of the wall moving at constant velocity is also 
to eliminate the less fitted individuals. Thus the size of 
the population can vary and even vanish. Besides, the 
models of evolution studied in [llj [2(| have a broader 
interest in physics as they are also related to models 
that appear in the theory of disordered systems [251 ] . in 
reaction-diffusion problems [26l . [27j or in particle physics 

Hi. mm. 

In the first part of this letter, we consider a continuous 
branching random walk in presence of a wall moving at 
constant velocity v (for the simulations below we will use 
a discrete space-time version to facilitate the program- 
ming). The diffusion is described as follows : during dt, 
the position of an individual is shifted by an amount r\\fdi 
where r\ is a random variable of zero mean and of variance 
2. For an initial individual at distance x > from the 
wall at t = 0, one defines the survival probability Q(x, t) 
that at least one of its descendants is still living at t (i.e. 
1 — Q(x,t) is the probability that all the descendance of 
the individual initially at position x ahead of the wall 
has been absorbed by the wall). To write the evolution 
of the probability Q(x,t), it is convenient to decompose 
the time interval [0, t + dt] into the first time interval dt 
and a remaining time interval [dt, t + dt]. During the first 
interval dt, while the wall moves by vdt, the single par- 
ticle initially at position x diffuses and can branch with 
probability (3dt into two particles, as expressed by the 
two terms of the following equation : 



Q(x + vdt,t + dt) = 



— — Q(x + riv dt, t) 

47T (1) 



+ 0dt(Q(x,t) -Q(x,t) 2 ) 
which becomes by taking the limit dt —> 
d t Q = d 2 x Q - vd x Q + (3(Q - Q 2 ) 



(2) 



with the boundary conditions Q(Q,t) — and Q(x,0) = 
1. It is important to notice that @ is exact. The origin 
of the term Q(x,t) 2 is that after a branching event, the 
two offspring have independent evolutions and therefore 
their survival probabilities are uncorrelated. 

All the x and t dependence of Q(x,t) can be extracted 
from the analysis of @ which is, up to the boundary 
conditions, the Fisher-KPP equation [ll|,[l2| m a moving 
frame at velocity v. 

Q{x,t) is obviously a monotonic function (increasing 
with x and decreasing with t : it increases if one starts 
further from the wall and it can only decrease with time). 
Therefore its t —* oo limit Q*(x) exists and satisfies 



dlQ* - vd x Q* + (3(Q* - Q* 







(3) 



which is the usual equation for the shape of a front solu- 
tion of the F-KPP equation moving at velocity v on the 
infinite line [lj|. 

On the infinite line [l6l |. for every velocity v, there is 
(up to a translation) a solution Q v (x) of (J3J) such that 
Q v (x) —> as x — > — oo and Q v (x) = I as x — > oo. 



These solutions are monotonic functions of x for v > v c = 
2->/p whereas for v < v c they have damped oscillations 
Q v (x) ~ Ae% x cos[-\/v% — v 2 (x — £o)/2] as x — > — oo. The 
relation between the decay of Q v (x) as x — ► — oo and the 
velocity is easy to obtain: one assumes that Q v (x) ~ e lx 
in the region where Q v (x) <C I and one gets from (|3|) 
that v is related to 7 by 



"(7) = 7 + 



7 ' 



(4) 



The monotonic decay for v > v c or the damped oscilla- 
tions v < v c are simply due to the fact that the solutions 
7 of the equation 1^(7) = v are real or complex conju- 
gate 7 = 7ft ± i'ji. It is also known that, at v = v c , the 
asymptotic decay has the form 



Q Vc (x) ~ (-A(x - x ) + B)e^ x - X ^ 



(5) 



where xq is arbitrary and 7 C = yp is the value of 7 for 
which 17(7) is minimum (and where v("f c ) = v c ) ■ 

With the boundary condition Q*(0) = the solution 
of is obviously Q*(x) — for v > v c as the wall moves 
faster than the spreading velocity v c of the population. 
On the other hand, for v < v c , the solution is the right- 
most positive arch of the solution on the infinite line. As 
v approaches v c , the region where Q* is small and where 
the quadratic term in can be neglected becomes larger 
and larger. For < v c — v <C 1 the solution Q* can 
be decomposed into two parts: a region x > L where 
Q*(x) is of order 1 and resembles the solution Q Vo (x) 
on the infinite line for v — v c (whose asymptotics decay 
is <j5j) ) and a region < x < L where Q*(x) is small 
and of the form Q*(x) — C sm(jix)e' lRX . To match the 
asymptotic decay of ([5]) near x = L, one should take 
L = 71-/7/ = 2-K/y/vl - v 2 and C = Ae~' fRL /j I such 
that 



*(x) ~ ALsin ( — 
V J_j 



(G) 



This implies that as v 



the survival probability 



vanishes with an essential singularity 



" (x) ~ exp 



7r V7c 



(7) 



One can remark that a shape of Q*(x) made up of two 
parts is almost exactly the same as the shape of a moving 
front on the infinite line, in presence of a cut-off [l7l ]. 

To obtain the long time dependence of Q(x,t) we as- 
sume that one can still decompose Q(x, t) into two parts : 
a region < x < L t where Q(x, i) is small with a solution 
of the form 



Q(x, t) = AL t sin 



. /7ra 
an — 



(8) 



(this assumption will be tested in the figure [3] below) and 
a region x > L t where Q(x, t) is not small and resembles 
([S]). One can check that (JSJ) is indeed solution of in 
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Velocity 


V < We 


V = V c 


V > Wc 


Q*(x) 


~exp( 7cA / 








Relaxation 




^ e -7c(37r 2 «"( 7c )/2) 1/3 i 1 /3 




Times r, f 









TABLE I: Critical behaviour (|7I12I13I15|I near v c of the survival probability Q{x,t) at large time t for an arbitrary v(j). 



the range < x < L t to leading non-vanishing order in 
L^ 1 and u — w c provided that L t satisfies 



Both (|T3j) and agree with earlier results by Kesten 
and Harris and Harris @ . 



d t L t = v 



(9) 



As one expects (at least for v < v c ) that L t /L ~ Lt(^c — 
v) 1 / 2 — > as i — > 0, this determines Lt and one gets that 
for < v e — v < 1, 



V7c(w c -v) V ^ 

where the function F(z) satisfies d z F 
solution (given that F(0) = 0) is 



vf'hj (10) 

-1 + F~ 2 whose 



- 2F + In 



l-F 



= 2z i.e. F = tanh(F 



(11) 



For large z, one gets F ~ 1 — 2e~ 2z ~ 2 . Therefore |TT 
implies that for w — > i;^ the convergence time r of the ex- 
ponential decay of Q(x, t) towards the fixed profile Q*(x) 



(v c - t-r 3/ V(2V7^) 



(12) 



For small z, one can expand the solution of (TTT]) 
F(z) = (Sz) 1 / 3 - 3z/5 + 0(z 5 / 3 ). This leads to L t ~ 
(37r 2 /7 c ) 1 / 3 t 1 / 3 + 0((u — u c )t) and therefore to a stretched 
exponential decay of Q(x, t) at v = v c 



Q(x,t)~eM-(^ 2 /j c ) 1/3 t 1/3 ] 



(13) 



At any finite t, one can notice (| 1011 1[) that L t can be 
expanded in powers of v — v c and one can get the whole 
time dependence for v > v c by an analytic continuation 
of the case v < v c . In particular, the scaling function F 
which replaces (|10l) for v — v c > is related analytically 
to the function F obtained for v — v c < and (jlOll 1[) are 
replaced by : 



L t =^L=F(^ {v - Vc 
F(z) =tan(F(z) - z). 



(14) 



Clearly F(z) = z + ir/2 - 1/z + 0(1/ z 2 ) for large z. 
Thus for < v - v c <C 1 and (v - v c ) 3 ^ 2 t > 1, one has 
L t ~ (v — v c )t and for fixed x : 



Q(x,t) ~ exp[7 c (x - (v - v c )t)}. 



(15) 



To study numerically the time evolution of the survival 
probability Q(x, t), it is more convenient to use a discrete 
space time version of the problem. Here we consider the 
case where x takes only integer values. At every time 
step, t — > t + 1, each individual has k offspring and each 
offspring has a probability p y of being produced at po- 
sition x + y. Then the time evolution of Q(x, t) which 
generalizes @ is given by 



Q(x + v,t+l) = l- 1 - ^2p y Q(x + y,t) 



(16) 



and by looking at the region where Q(x,t) ~ e lx is small, 
the expression ((4| of the velocity v(-f) becomes 



v(i) = - fog 

7 



~y 



As before, the spreading velocity v c of the population in 
absence of the wall is v c = min 7 ^(7) . In this case one can 
generalize the above calculations to an arbitrary u( 7 ) and 
the results (|TI 1 21 1 31 1 5[) are replaced by the expression of 
tabled] 

To facilitate the analysis and avoid periodic or ape- 
riodic dependence due to the commensurability of v we 
decided to fix v = 1 for the velocity of the wall and to vary 
v c by changing the p y . In all the figures below, we show 
the case where P2 = 1/18 + 6,p\ = 1/18, po = 16/18 — 5. 
Then by varying 5 we can vary v c and therefore the dif- 
ference v c — v (which vanishes for 5 = 0). More pre- 
cisely, for small v c — v, one has the first-order relation 
6 ~ I5(v c — w)/(21n4) in the present case. 

In figure [3] we test our assumption ([5]) . For an initial 
individual at x — 1, i.e. for Q(x, 0) = 5 x i, we computed 
Q(x,t) by iterating JTgJ). We plot Q(x,t)e" l{x - L ^ j L t 
which according to ([S]) should depend on x/L t only. In 
figure [H we compare the values of Q(l,t) obtained by 
iterating p^]) with the theoretical predictions (JSJ [TO]). 
The agreement is excellent. 

At t — 00 , the size of the population is cither zero 
(with probability 1 if v > v c and 1 — Q(x, 00) if v < v c ) 
or infinite. At large but finite final time T, there are 
however events corresponding, for example, to having 
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FIG. 2: The points represent the logarithm of the survival probability Q(x, t) at x = 1 for different values of v — v c both below 
and above the critical velocity, obtained by iterating (I16|l . The lines are the theoretical predictions (|8I10[1 . 




FIG. 3: Shape of the rescaled survival probability 
Q(x,t)e (c ^ x ~ L *> /Lt in the linear domain at different times 
t = 10 5 (dashed), 5.10 5 (circle), 11.10 s (line). The cor- 
responding lengths are L t = 1.05.10 2 , 2.01.10 2 , 2.93.10 2 ) for 
v — v c ~ 10~ 4 > (5 = 2~ntl( Vc ~ v ) at nrst or< ier for this 
particular model). The length L t is measured directly as the 
point where Q(Lt,t) = 0.5 by linear interpolation between 
the sites of the lattice. The observed shape agrees with (JSJ. 



exactly one survivor. One can then try to understand 
which strategy leads to such events. For v < v c , t 1 
and T—t = t' 3> 1, the system adopts a quasi-stationary 
regime during which its properties do not depend on the 
initial conditions (number and positions of the individu- 
als at t = 0). In particular, the average population size 
at time t conditioned by the survival of only one individ- 
ual at final time T = t + 1', do not depend on t as long 
as t > 1 and t' > 1. 

We are now going to sketch how the properties of this 
quasistationnary regime can be computed from ^ or 
(TIT)]) . Details will be published in a future work [31| ■ Let 



us introduce the generating function 



G(x,t;f) = (l[e 



,-/(*") 



(17) 



where / is an arbitrary function, Nt is the size of the 



population at time t and x\ 



„(*) 

C N t 



are the positions of 



individuals at time t. By analysing what happens during 
the very first instant dr, as in (fT]), one finds that : 

G{x + vdt, t + dt; f) =f3dt(G{x, t; f) 2 - G(x, t; /)) 

e-^^dri 



-in 



G(x + r]Vdt,t;f). 

(18) 



Therefore 1 - G(x, t; f) satisfies the F-KPP equation (JSJ) 
with the boundary conditions G{x — 0,t\ /) = 1 and 
G(x, t = 0; /) = e~f( x \ In the discrete space time case, 
a similar reasoning shows that 1 — G satisfies (Tl7jl) . 

One can also consider generating functions at two (or 
more) different times < t < t + 1' defined as : 



t+t' 



-h(x 



(t+t 1 ). 



4=1 j = l 

(19) 

Like in (p~5|) . one can show that 1 — H (x, t, t'; /i, f^) 
as a function of x and t satisfies (fT]) or (TT51) with 
boundary conditions H(x = 0, t, t'; /i, /2) = 1 and 



H{x,0,t';hJ 2 ) 



-/i(x) 



G{x,t' ; $2) which, as we have 



just seen, can also be determined from (TTJ) or (TTo|) . 

If one chooses /1 = A and ji — /i, then 
H(x, t, t'; fx, f 2 ) = ( e - XNt ~^ N W) is the generating func- 
tion of the sizes of the population at two different times. 
If one expands H to first order in powers of e _M and of A 



H(x, t, t'; A, fx) ~H 00 (x, t, t') + e~^H w {x, t, t')+ 
XH 01 (x, t,t') + (x, t, t') 



(20) 
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one gets that the average size of the population 
(N t }x=t+f at time t conditionncd on having a single sur- 
vivor left at time T = t + t' is given by : 

{N t ) T =t+t> = -ff u (i, t, t')/H w (x, t, t'). (21) 

Numerical values obtained by iterating the equations 
satisfied by Hqq, -Hio, Hqi and H\\ (derived by replacing 
H by its expansion ([2U|) into (HE])) are shown in figure 0$. 
By pushing further the expansion (|20[) in powers of A, one 
can have access to higher moments of the population size 
N t . 

In figure [H we see that for v < v c both (N t }T and 
(N£)t I '{NtlT become constant over a very long time in- 
terval, indicating that there is a quasi-stationary regime. 
As d — > v c , the size of the population in this quasi- 
stationary regime diverges and the ratio {Nf)T/{N t )\ 
seems to converge to a value close to 2. This indicates 
that the fluctuations of the population size Nt in the 
quasi-stationary regime are of the same order as (Nt) 
itself. 

For v < v c in fact, the whole distribution of Nt at 
<C t conditionned on the survival of a single individual 
at T = t + t 1 > t is independent of the position x of the 
initial individual and reaches a quasi-stationary regime 
@ before a final relaxation to the final state containing 
only one individual. In this quasi-stationary regime one 
can show from and (|20I21[) that the average size of the 
population (N t ) qs can be calculated for v — > v c (details 
will be published in [3l|) and one finds 



(N, 



t/qs 



exp 



ir 2 v"(-f c ) 
2(v c - vj 



(22) 



For v > v c , no quasi-stationary state is observed (see 
figure . Instead, the evolution of the system condi- 
tionned on having a single survivor at a later time T can 
be divided into two parts. The first one corresponds to a 
rapid growth of the population, which can be described 
by the modified dynamics introduced in [f|, followed by 
a regular absorption of individuals by the wall leading to 
the survival of exactly one individual at T. Numerical 
results obtained from (fTo| and (|20I21[) are presented in 



figure [5]: they indicate that (N t )r takes a scaling form 
(N t ) T ~ exp[T/(t/T)]. 

In this article, we have seen how a branching random 
walk in presence of a moving absorbing walk undergoes 
a phase transition at a critical velocity v c with a survival 
probability which vanishes with an essential singularity 
([7]) and decays with time like a stretched exponential (fT3"|) 
at the transition. We have also seen how the properties of 
the quasi-stationary regime § can be computed (fTTH2"Tj) 
from the travelling wave equation. 

It is interesting to notice that several of our results are 
very reminiscent of what is known of noisy F-KPP like 
equation [13, . First the relation (f2"2")l between the size 
TV of the population in the quasi-stationary regime and 
the velocity can be rewritten as : 



7rV'( 7c ) 1 
2 7c 2 In 2 N 



(23) 



which is exactly the prediction of the cut-off approxima- 
tion [It], US, Moreover, for v < v c , the relaxation 
time r ~ (v c — u)~ 3 / 2 can be rewritten as : 



ln 3 iV 



(24) 



which is exactly what was observed recently for the co- 
alescence times in the genealogies of evolution models 
[2 (J . |33j related to these noisy travelling wave equations. 

Our results on the quasi-stationary regime presented 
here (eq. (|22p and figures 2] and [5]) are preliminary and 
will be developed in a forthcoming work [3lj . 

Beyond the understanding of this quasi-stationary 
regime, it would be interesting to see whether our main 
results which are based on our guess ([5]) for the shape of 
the solution of © could be recovered by more standard 
methods like the renormalization group approach Q or 
the phenomenological picture [HJ of noisy fronts. Among 
other possible extensions of the present work, one could 
try to study how our results would be modified when the 
position of the wall has a more complicated evolution 
than a constant velocity [34| . 
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FIG. 4: On the left side : average size of the population at < t < T conditionned by the existence of a unique survivor 
at T — 250000 for different values of v — v c < 0. The results are obtained using fj 16120 12 1 11 . On the right side, the ratio 
(N?} T /(Nt}T is represented for < t < T for the same values of v — v c < 0. In both cases, the initial individual at t = is 
at x = 1. A quasi-stationary regime (with constant (N t )T and {N%)t / {Nt)'r) is observed as soon as T is much larger than the 
relaxation time r. 
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FIG. 5: The logarithm of the average size of the population 
divided by T at < t < T conditionned by the existence of a 
unique survivor at T is plotted as a function of the rescaled 
time t/T for v — v c = 10~ 2 and different final times T — 
62500, 125000, 250000. Unlike figure H no quasi-stationary 
regime is observed. 
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